Elliptic flow in Pb+Pb collisions at y/sj^ = 2.76 TeV: 
hybrid model assessment of the first data 
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We analyze the elliptic flow parameter V2 in Pb+Pb collisions at y^sjvjv ~ 2.76 TeV and in 
Au+Au collisions at ^Jsnn = 200 GeV using a hybrid model in which the evolution of the quark 
gluon plasma is described by ideal hydrodynamics with a state-of-the-art lattice QCD equation of 
state, and the subsequent hadronic stage by a hadron cascade model. For initial conditions, we 
employ Monte-Carlo versions of the Glauber and the Kharzeev-Levin-Nardi models and compare 
results with each other. We demonstrate that the differential elliptic flow V2{pt) hardly changes 
when the collision energy increases, whereas the integrated V2 increases due to the enhancement of 
mean transverse momentum. The amount of increase of both V2 and mean pr depends significantly 
on the model of initialization. 
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The recently started heavy ion program at Large 
Hadron Cohider (LHC) in CERN opens up opportuni- 
ties to explore the deconfined matter, the quark gluon 
plasma (QGP), in a wider temperature region. Elliptic 
flow which played an essential role to establish the 
new paradigm of the strongly coupled QGP 0, iSj] at Rel- 
ativistic Heavy Ion Collider (RHIC) in Brookhaven Na- 
tional Laboratory (BNL) is one of the key observables 
at LHC to investigate the bulk and transport properties 
of the QGP. First elliptic flow data in Pb-|-Pb collisions at 
^SNN = 2.76 TeV were recently published by the ALICE 
Collaboration Q . The first goal of flow measurements is 
to see whether hydrodynamic models reproduce the flow 
as well at LHC as at RHIC, and thus whether the QGP 
depicts similar strong coupling nature at LHC. 

This Rapid Communication is a sequel to our previous 
work where we predicted the elliptic flow parameter 
1)2 before any LHC data was available. In this publica- 
tion we take the advantage of the first LHC data [H, 
to fix the final particle multiplicity, which removes the 
main uncertainty in our prediction, and allows us to use 
a Glauber type initialization too. We calculate the el- 
liptic flow parameter V2 and its transverse momentum 
{pt) dependence in Pb-|-Pb collisions at LHC and com- 
pare them with the data. Our model for the space-time 
evolution of the matter is the same we used in Ref. Q: 
A hybrid model where the expansion of the QGP is de- 
scribed by ideal hydrodynamics Q, and the subsequent 
evolution of hadronic matter below switching tempera- 
ture Tsw = 155 MeV, is described using a hadronic cas- 
cade model JAM 0. During the fluid dynamical stage, 
we employ EoS s95p-vl.l, which interpolates between 



hadron resonance gas at low temperatures and recent lat- 
tice QCD results by the hotQCD collaboration [Tol . [Til 
at high temperatures in the same way as s95p-vl [12l |. 
but the hadron resonance gas part contains the same 
hadrons and resonances as the JAM hadron cascade 0]. 
The details of the interpolating procedure arc explained 
in Ref. [12] and the parametrization and EoS tables are 
available at Ref. [T3| . 

Initial time of hydrodynamic simulations is fixed to be 
To = 0.6 fm/c throughout this work. For initial condi- 
tions in the longitudinal direction, we assume the Bjorken 
scaling solution To initialize the density distribu- 

tions in the transverse plane, we utilize two Monte-Carlo 
approaches: Monte-Carlo Glauber (MC-Glauber) model 
[la l and Monte-Carlo Kharzeev-Levin-Nardi (MC-KLN) 
model [l^. Using these Monte-Carlo models, we cal- 
culated initial conditions for hydrodynamic simulations 
in the transverse plane with respect to the participant 
plane in our previous work . These initial density pro- 
files contain effects of eccentricity fluctuation on average. 
However, the ALICE Collaboration mainly obtained V2 
using the 4-particle cumulant method W2{4} [13 1 which 
is expected to contain less eccentricity fluctuation and 
non-flow effects than V2 measured using the 2-particle 
cumulant method, W2{2}[3- Therefore we calculate in 
this Rapid Communication initial profiles with respect to 
the reaction plane: We average over many events using 
Monte-Carlo calculations instead of shifting and rotating 
a distribution event- by- event to match the main and sub 
axes of the ellipsoids as was done in the previous work 
[^ [131 ■ It should be noted that the distributions obtained 
in this way are not identical to the ones from the optical 
Glauber model or the factorized KLN (fKLN) model [U 
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due to finite nucleon size effects [l6[ : tfie collision points 
in the transverse plane are smeared using inelastic cross 
section oi p + p collisions in the "mean-field" option in 
the Monte-Carlo approach [l^ to obtain smooth initial 
conditions for hydrodynamic simulations. 

In the MC-KLN model, we calculate distribution of 
gluons at each transverse grid using the fct-factorizcd 
formula [l^. Using the thickness function T^, we 
parametrize the saturation scale for a nucleus A as 
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and similarly for a nucleus B. We choose A = 0.28 and a 
proportionality constant in the unintegrated gluon distri- 
bution in the fct-factorized formula to reproduce ccntral- 
ity dependence of pt spectra obtained by the PHENIX 
Collaboration As a default parameter set at LHC, 
we use the same parameters except for colliding energy 
and mass number of incident nuclei. This predicted 
dNch/drj ~ 1600 at 5% most central collisions [6|, which 
turns out to be consistent with the recent ALICE mea- 
surement 0, nil- 

In the MC-Glauber model, one calculates the number 
distributions of participants Ppart and of binary collisions 
PcoU for a given nuclear density distribution. We model 
the initial entropy distribution in hydrodynamic simu- 
lations as a linear combination of ppart and pcoU hi the 
transverse plane: 
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At the RHIC energy, the mixing parameter a = 0.18 
and the proportionality constant in Eq. ([2]) are chosen 
to reproduce the centrality dependence of pr spectra at 
RHIC [1^. We tune these two parameters in Pb-|-Pb 
collisions at LHC to reproduce the centrality dependence 
of charged hadron multiplicity [2^. For both initializa- 
tions we do the centrality cuts accordingto the A^part 
distribution from the MC-Glauber model Q 

In Fig.[TJ we calculate dNch/drj/ (A^part/2) as a function 
of A'part for initial conditions from the MC-Glauber and 
the MC-KLN models and compare them with data. The 
experimental data point in inelastic p + p collisions at 
y/IWN = 2.36 TeV ^ is plotted at A^part = 2. The MC- 
KLN initialization leads to remarkable agreement with 
the ALICE data. On the other hand, it is difficult to fit 
the data within the current two-component picture in the 
MC-Glauber model: The results from the MC-Glauber 
initialization with a = 0.08 almost trace the ones from 
the MC-KLN initialization and the ALICE data except 
for 0-5% and 70-80% centrality 

Shown in Fig. [2] is the initial eccentricity with respect 
to reaction plane as a function of Apart in Pb-f Pb col- 
lisions at y^SNN = 2.76 TeV and in Au-|-Au collisions 
at y^SNN = 200 GeV. As previously known, the kt- 
factorized formula of KLN model generates larger ec- 
centricity than the Glauber model does [U [2^. In 
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FIG. 1. (Color online) Centrality dependence of charged 
hadron multiplicity in the MC-Glauber and the MC-KLN ini- 
tialization is compared with ALICE data [S^, H^. A data 
point from inelastic events at ^sjvjv ~ 2.36 TeV in p + p 
collisions [2^ is shown at A^part = 2. Each point in theoreti- 
cal results from right to left corresponds to 0-5, 5-10, 10-20, 
20-30, 30-40, 40-50, 50-60, 60-70, and 70-80% centrality, re- 
spectively. 
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FIG. 2. (Color online) Eccentricity with respect to the re- 
action plane as a function of A^part in Pb-f Pb collisions at 
^ysNN = 2.76 TeV and in Au-|-Au collisions at y^ijvjv ~ 200 
GeV. Each point from right to left corresponds to 0-5, 5-10, 
10-20, 20-30, 30-40, 40-50, 50-60, 60-70, and 70-80% central- 
ity, respectively. 



the MC-KLN model, eccentricity in Pb-|-Pb collisions 
at y/SNN = 2.76 TeV is slightly larger than that in 
Au-I-Au collisions at ^/snn = 200 GeV when the central- 
ity is fixed 0]. On the other hand, in the MC-Glauber 
model, eccentricity in Pb+Pb collisions at y^SNN = 2.76 
TeV is slightly smaller than that in Au-t-Au collisions at 
y^SNN = 200 GeV for a fixed centrality. 

This is due to the smearing process we use to obtain 
a smooth initial profile for hydrodynamic evolution. As 
mentioned, we use the inelastic cross section in p + p col- 
lisions, CTin, to smear the distribution of collision points. 
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FIG. 3. (Color online) Transverse momentum distribution of charged hadrons at 10-20% (circles) and 40-50% (squares) 
centralities in Pb-|-Pb collisions at ^/sNN = 2.76 TeV (filled symbols) and in Au+Au collisions at ^/snn = 200 GeV (open 
symbols). Results from (a) the MC-Glauber initialization and (b) the MC-KLN initialization. For the sake of comparison and 
visibility, the spectra are scaled by 2, 1/10, and 1/5 for 10-20% at RHIC, 40-50% at LHC, and 40-50% at RHIC, respectively. 



This cross section is ^ 1.5 times larger at LHC than 
at RHIC, and thus the smearing area, 5* = a-m [l6| . is 
also larger at LHC, and the eccentricity is reduced. Our 
smearing procedure also leads to a smaller eccentricity 
than the conventional value of MC-Glauber mode0. The 
effect of smearing is smaller in the MC-KLN initializa- 
tion, and we have checked that the eccentricity at LHC 
turns out to be essentially the same as at RHIC when 
the smearing area is the same. Systematic studies of 
initialization and its effects will be shown in a later pub- 
lication [23. 

Figure [3] shows comparison of transverse momentum 
distributions of charged hadrons between RHIC and LHC 
energies at 10-20% and 40-50% centralities. As clearly 
seen from figures, the slope of pt spectra becomes flatter 
as collision energy and, consequently, pressure of pro- 
duced matter increases. To quantify this, we calculate 
mean pT of charged hadrons. In the MC-Glauber initial- 
ization, mean pT increases from RHIC to LHC by 21% 
and 19% in 10-20% and 40-50% centrality, respectively. 
On the other hand, the corresponding fractions are 25% 
and 24% in the MC-KLN initialization. Since our calcu- 
lations at RHIC were tuned to reproduce the p^-spectra, 
this means that at LHC the spectra calculated using the 
MC-KLN initialization are slightly flatter than those cal- 
culated using the MC-Glauber initialization. 

We compare integrated V2 for charged hadrons with 
resp ect to reaction plane with the ALICE Q and STAR 
[28| W2{4} data in Fig.|4l When evaluating the integrated 
V2, we take account of both transverse momentum and 
pseudorapidity acceptance as done in the experiments, 
i.e. 0.2 <pT < 5.0 GeV/c and | 7? |< 0.8 for ALICE, and 



0.15 <PT < 2.0 GcV/c and | r; |< 1.3 for STAR. We want 
to emphasize that not only the pT cut [2^ , but also the 
pseudorapidity cut plays an important role in a consistent 
comparison with the data. Due to the Jacobian for the 
change of variables from rapidity y to pseudorapidity 77, 
V2{y = 0) < V2{ri = 0) for positive elliptic flow 30lpl. In 
the case of the MC-Glauber (MC-KLN) initialization in 
40-50% centrality, V2 integrated over the whole pT region 
is -14% (-10%) larger at 77 = than at y = 0. 

When the MC-Glauber model is employed for initial 
profiles, centrality dependence of integrated V2 from the 
hybrid approach almost agrees with both ALICE and 
STAR data. Since eccentricity fluctuation contributes lit- 
tle and negatively to W2{4} in non-Gaussian distribution 
of eccentricity fluctuation [l^, [l^, this indicates there 
is only little room for the QGP viscosity in the model 
calculation. On the other hand, apparent discrepancy 
between the results from the MC-KLN initialization and 
the ALICE and STAR data means that viscous correc- 
tions during the hydrodynamic evolution are required. 

From RHIC to LHC, the pr-integrated V2{\ V \< 0-8) 
increases by 24% and 25% in 10-20% and 40-50% central- 
ity, respectively, in the MC-Glauber initialization. On 
the other hand, in the MC-KLN initialization, the in- 
crease reaches 42% and 44% in 10-20% and 40-50% cen- 
trality, respectively. Since eccentricity does not change 
significantly (at most ±6% in 40-50% centrality) from 
RHIC to LHC as shown in Fig. [51 the significant increase 
of integrated V2 must be attributed to a change in trans- 
verse dynamics. 

Finally, we compare 7^2 (p t) of charged hadrons with 
ALICE H and STAR data in 10-20% (Fig. [D^a)) 



In the MC-Glauber model in the literature Il5l . one assumes 
5 function profile for each collision point in Ppart distribution 
rather than a box-like profile in the present work. 



^ Notice that even if one assumes the Bjorken scaling solution, one 
has to consider the pseudorapidity acceptance since V2{ri) is not 
constant even if V2 {y) is [30| . 
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FIG. 4. (Color online) Centrality dependences of vi for charged hadrons with respect to reaction plane (a) in Pb+Pb collisions 
at ^/s^ = 2.76 TeV (| |< 0.8, 0.2 < pr < 5 GeV/c) and (b)_in Au+Au collisions at ./t^ = 200 GeV (| t] |< 1.3, 
0.15 < pt < 2 GeV/c) are compared with ALICE and STAR [2^ v-2 data, respectively. ALICE data points are shifted 
horizontally for visibility. 
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FIG. 5. (Color online) Transverse momentum dependences of V2 for charged hadrons in the MC-Glauber (circles) and the 
MC-KLN (squares) initialization are compared with ALICE [1| (triangles) and STAR ^ (band) W2{4} data in (a) 10-20% 
centrality and (b) 40-50% centrality. We take account of pseudorapidity cut, | r\ |< 0.8 (1.3), in the ALICE (STAR) data. 



and 40-50% (Fig. [5jb)) centrality. Interestingly, the data 
at LHC agrees with the data at RHIC within errors. The 
calculated vi {pr ) shows similar independence of collision 
energy when MC-Glaubcr initialization is used, whereas 
MC-KLN initialization leads to a slightly larger V2{pt) 
at the larger energy. For MC-Glauber results, the fit to 
data is fair below pT ^ 1.5 GeV/c and pt ~ 0.8 GeV/c 
momenta in the 10-20% and 40-50% centralities, respec- 
tively. Results from the MC-KLN initialization at both 
energies are significantly larger than experimental data 
in the whole pt region, which again indicates necessity of 
viscous corrections in hydrodynamic evolution. For both 
initializations the difference between the data and the 
calculated V2{pt) is larger in more peripheral collisions. 
This too can be understood as an indication of viscosity, 
since the more peripheral the collision, the smaller the 
system and the more anisotropic its shape, and both of 
these qualities enhance the dissipative effects. 



Due to the relationships among the px spectrum, px 
averaged V2, and px differential V2{px), the flatter the px 
spectrum, the larger the V2 even if 1)2 (pt) stays the same. 
It is also worth noticing that the steeper the slope of 
V2{px), the larger the increase in V2 for the same increase 
in mean px ■ This is the main reason why quite a similar 
increase of mean px for both MC-Glauber and MC-KLN 
initializations leads to much larger increase of V2 for MC- 
KLN than for MC-Glauber initialization. 

At the time of this writing, the initial state of the fluid 
dynamical expansion of heavy-ion collisions at ultrarela- 
tivistic energies is quite uncertain. This has been a long- 
standing issue in the physics of heavy ion collisions which 
must be by all means resolved. If color glass condensate 
(CGC) [3l| initial conditions, like the ones obtained us- 
ing the MC-KLN model in the present work, are realized 
in nature at both RHIC and LHC energies, the larger 
deviation of V2 from the data at LHC than at RHIC in 
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Figs, m and [5] could mean that viscous effects are larger 
at LHC than at RHIC. This can indicate a larger spe- 
cific shear viscosity, rj/s, at larger temperatures. For a 
better interpretation of current experimental data, the 
issue should be clarified in near future by determining 
the initial conditions better and by a more detailed anal- 
ysis using a hybrid model of viscous hydrodynamics and 
hadron cascade [s^ . 

In summary, we calculated transverse momentum dis- 
tribution of charged hadrons, centrality dependence of 
integrated elliptic flow parameter V2, and differential el- 
liptic flow V2{pt) in Pb+Pb collisions at y^s^N = 2.76 
TcV and in Au+Au collisions at ys^vjv = 200 GcV. We 
compared V2 and V2{pt) with respect to reaction plane 
from the hybrid model with V2 data mainly obtained 
from the 4-particle cumulant method. Transverse mo- 
mentum distributions become harder, whereas the shape 
of V2{pt) does not change so much as the collision en- 
ergy increases. Thus the increase in pr-intcgratcd V2 is 
due to the increase in mean pT- However, the intrinsic 
slope of V2{pt) depends on the initialization: The slope 
from the MC-KLN initialization is steeper than that from 
the MC-Glauber initialization, and thus essentially the 
same change of mean pT leads to larger increase of pr- 
integrated V2 for MC-KLN initialization than for MC- 



Glauber initialization. The larger difference between the 
data and our MC-KLN result at LHC than at RHIC may 
indicate larger dissipative effects at LHC than at RHIC. 
All this emphasizes the importance of understanding ini- 
tial conditions in relativistic heavy ion collisions towards 
extracting the bulk and transport properties of the QGP. 
In future, it would be interesting to compare our results 
with data obtained using a more sophisticated elliptic 
flow analysis [l9j . in which both non-flow and eccentric- 
ity fluctuation effects are removed. 
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